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We investigate the consequences of the stripe formation in the nickelate La5/3Sri/3Ni04 for the details of 
its crystal structure and electronic states. Our data are based on numerical simulations within density func- 
tional theory (DFT) and the generalized gradient approximation (GGA). The on-site Coulomb interaction is 
included in terms of the LDA-l-U scheme. Structure optimization of preliminary experimental data indicates 
a strong interaction between the structural and electronic degrees of freedom. In particular, we find a segre- 
gation of the diagonal filled stripes induced by a delicate interplay with the magnetic coupling. Beyond the 
cooperative effect of stripe segregation and spin order, distinct octahedral distortions are essential for the 
formation of an insulating state. 



1 Introduction 

Doped Mott insulators present a wealth of intriguing properties with a great potential for technological 
applications. They are either based on metalicity, like high-Tc superconductivity and thermoelectricity, or 
trace back to an insulating behavior, like colossal magneto-resistance and multiferroics. Therefore, un- 
derstanding the origin of the insulating state of these doped Mott insulators has clear implications, be it 
caused by a lack of coherence [1] or by a formation of an inhomogeneous state, as a stripe phase, for in- 
stance [2-8] . Though intensely debated, see Ref. [9] for a review, there is growing consensus that electronic 
models can account for the poorly metallic character of the diagonal filled stripes observed in cuprates, al- 
most irrespective of their microscopic structure. For instance, magnetic diagonal filled stripe phases of the 
Hubbard model are insulating [10], while a poorly metallic phase with nodal quasiparticles characterizes 
flux diagonal filled stripe phases of the t — J model [11]. In contrast, similar model calculations devoted 
to isostructural nickelates fail to reproduce an insulating groundstate [12]. 

Experimentally, the stripe order in La2_xSrxNi04 and Nd2_xSrxNi04 has been evidenced by spin- 
polarized as well as unpolarized neutron scattering. With decreasing temperature, the appearance of peaks 
at the wave vectors — 27r(e, e) and Q^. — 7r(l ±e, lie), respectively, with e ~ x, is interpreted as onset 
of a diagonal filled stripe. The onset temperature is sensitive to the doping, turning out to have a maximum 
at x ~ 1/3, i.e. when the charge wave vector matches right the spin wave vector In this case, specific 
heat [13], transport [14], and optical conductivity [15] data point to an ordering temperature of T ~ 240 K. 
Recently, the stability of the stripe phase has been further probed via non-linear transport experiments for 
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Fig. 1 Scliematic representation of tlie Ni02-layer in La5/3Sri/3Ni04 iiighiligiiting tlie diagonal Ni stripes. O atoms 
are not shown for clarity. Arrows indicate the principal structural distortion consistent with AF spin order. 

the system Ndi.67Sro.33Ni04 [16]. No depinning of the stripes has been observed: instead, the non-linear 
transport is attributed to resistive heating, thereby pointing towards a rather robust insulating state. 

On the theoretical side, quite some efforts have been made to properly describe these materials, espe- 
cially by the investigation of realistic models. While an extensive study of a two-band Hubbard model for 
the Sg electrons correctly predicts diagonal filled stripes to be lowest in energy within this class of states, 
the groundstate is still metallic [12]. Adding electron-lattice coupling promotes the insulating groundstate, 
but the available results are limited to the x — 1/2 case and to particular lattice distortions [17]. Band 
structure calculations also have been performed [18, 19]. However, because they make use of an explicit 
lattice structure as input, which so far could not be determined experimentally for the title compound, the 
validity of the calculations could be limited. Indeed, as argued in [17], a strong spin-charge-lattice cou- 
pling is to be anticipated in layered nickelates. Using an two-band Hubbard model, Hotta and Dagotto [20] 
have demonstrated the interplay of Coulombic and phononic interactions and its relevance for the stripe 
formation. 

Another appealing scenario of spin-lattice interaction is sketched in Fig.[Tl where we show the Ni atoms 
for a Ni02 -layer Atoms along the stripes and in the magnetic domains are represented by full and open 
circles, respectively. A Ferromagnetic (F) as well as an antifeiTomagnetic (AF) order between the spins on 
the stripe and the domain sites is compatible with the experimental modulation vector = 7r(2/3, 2/3). 
For the AF alignment, the structural distortion shown in Fig.[T]enhances the superexchange over the stripes 
significantly, while it reduces the frustrated superexchange within the magnetic domains only moderately. 
The configuration thus is expected to be subject to an instability against a modulation of the Ni-Ni bond 
length perpendicular to the stripe direction. In contrast, the F alignment in this respect is consistent with 
an equal spacing of the Ni atoms, because a modulation of the bond lengths tends to reduce the kinetic 
energy. The competition between the two states makes it highly desirable to understand in detail the shape 
of the lattice deformations encountered in the system, and their consequences for the electronic structure. 



2 Technical details 

To that aim we have performed band structure calculations for the a; = 1/3 compound La5/3Sri/3Ni04. 
Our following considerations first focus on the lattice deformations tracing back to the stripe formation, 
where the parent body-centered tetragonal K2NiF4 structure with space group lA/mmni will act as a 
reference frame [21]. Lattice deformations have not been reported in the literature so far, because the 
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abbreviation 


intraplane spin alignment 


interplane spin alignment 


ADFAS 


F 


AF 


CDFAS 


AF 


F 


GDFAS 


AF 


AF 



Table 1 Definitions of the A-, C-, and G-type D(iagonal) F(illed) A(ntiferromagnetic) S(tripe) spin configurations. 

experimental refinement of the crystal structure has concentrated on the lattice parameters. In fact, refine- 
ment of the atomic positions is difficult without a proper structure model, which, finally, is provided by our 
structure optimization data. As second step, we analyze the electronic properties of the system for differ- 
ent spin configurations, which allows us to determine the magnetic groundstate of various diagonal filled 
antiferromagnetic stripes. The spin configurations under investigation are defined in table [T] In addition, 
the CDFAS and GDFAS patterns are illustrated in Fig.|2] 

In our electronic structure calculations we use density functional theory (DFT) as implemented in the 
Wien2k code [22], The latter is appropriate for the present investigation, because it allows us to perform a 
full structure optimization for a complex material with large unit cell [23] and to account for the electronic 
correlations present in the system, by the LDA+U (local density approximation plus on-site Coulomb 
interaction) method [24]. For the structure optimizations the generalized gradient approximation (GGA) 
[25] and for the LDAh-U calculations the self-interaction correction (SIC) [26, 27] is used. 

The tetragonal La2-xSrxNi04 unit cell resulting from the parent lA/mmm lattice contains two for- 
mula units. It is very convenient to describe the atomic arrangement in terms of parallel layers of NiOg- 
octahedra, separated by the La/Sr atoms. To be specific, we have Ni02-layers parallel to the a6-plane of 
the tetragonal unit cell, consisting of a square planar lattice of Ni atoms, and O atoms located midway 
between them. Due to the stripe charge order, we have two cry stallographic ally inequivalent Ni02-layers 
comprising the metal sites Nil/Ni3 and Ni2/Ni4, respectively. For further details on the crystal structure, 
the supercell setup, as well as other technical aspects, we refer the reader to Ref. [19]. We use the setup 
described therein as the starting point of our present study. However, our new calculations thoroughly take 
into account the structure relaxation induced by the stripe charge order. To this end, we optimize the atomic 
positions individually for all spin patterns under consideration. 




Fig. 2 Magnetic structure of the Ni atoms in the Ni02-layers of La5/3Sri/3Ni04 for the CDFAS (left hand side) and 
GDFAS (right hand side) spin configurations. Full yellow (blue) circles denote the Nil (Ni3) atoms, whereas dashed 
circles refer to the second Ni02 -layer, containing the Ni2 (Ni4) atoms, see the text. 
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3 Results and discussion 

Turning to the data obtained from our structure optimizations, we first mention that the NiOe-octahedra 
of La5/3Sri/3Ni04 are elongated and tilted with respect to the c-axis of the tetragonal unit cell [28,29]. 
While structure distortions, in general, are tiny for the ADFAS pattern, we observe a sizeable response of 
the NiOg-octahedra to the charge order for both the CDFAS and GDFAS pattern, i.e. for AF intraplane 
magnetic coupUng. The fact that the latter distortions are closely related to each other supports our initial 
speculations about a strong spin-lattice interaction. We summarize the observed apical and equatorial Ni-O 
bond lengths in the four inequivalent NiOg-octahedra in Table|2l revealing significant deviations from the 
ideal values of 2.21 A and 1.92 A, respectively. The distortion associated with these values is illustrated in 
Fig. |3]for the Ni2/Ni4-plane, where the shifts off the symmetric positions have been magnified by a factor 
of 4. The pattern is similar for the Nil/Ni3-plane, which is obtained by shifting the Ni2/Ni4-plane by the 
vector (a/2,0,c/2). While for stripe ordered Lai 875Bao.i25Cu04 [30], for example, a modulation of the La 
ions is reported, such displacements are negligible in our case for both the La and Sr sites. 

Whereas a strong relaxation of the structure has to be expected for a NiOg -octahedron, we additionally 
observe a characteristic relaxation pattern for the Ni sublattice, see the schematic illustration in Fig.[T] With 
respect to the (highlighted) diagonal stripes, neighboring Ni atoms are subject to displacements towards 



atom 


apical Ni-O bond length (A) 


equatorial Ni-O bond length (A) 


Nil 


2.17 (2x) 


1.94 (4x) 


Ni2 


2.23 (2x) 


1.94 (4x) 


Ni3 


2.20 (2x) 


1.89 (2x) and 1.92 (2x) 


Ni4 


2.24 (2 X) 


1.88 (2x) and 1.94 (2x) 



Table 2 Bond lengths in the NiOe-octahedra obtained for AF intraplane magnetic coupling (CDFAS/GDFAS). 




Fig. 3 Optimized crystal structure of a Ni02-layer in La5/3Sri/3Ni04. The Ni2 (Ni4) and O sites are represented by 
black (green) circles and diamonds, respectively, and stripes by dotted lines. 
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displacement (A) 



atom moving 


towards 


ADFAS 


CDFAS 


GDFAS 


Ni3 


Nil chain 


0.0002 


0.0065 


0.0016 


Ni4 


Ni2 chain 


0.0032 


0.0114 


0.0300 



Table 3 Atomic displacements as obtained from the structure optimizations for the different spin patterns. The values 
refer to the relative shift of a Ni3/Ni4 site towards the chain formed by the Nil/Ni2 sites within a Ni02 -layer. 



them. To be specific, all shifts are oriented perpendicular to the stripe direction and parallel to the basal 
plane of the tetragonal unit cell, thus parallel to the Ni02-layers. Their absolute values are summarized in 
Table[3] The distortion pattern of the Ni sublattice is well described in terms of a segregation of neighboring 
stripes. By this mechanism, the magnetic intraplane interaction is enhanced within the magnetic domains. 
After all, this is advantageous for the AF configuration only, since the F configuration is predominantly 
promoted by local processes. Furthermore, our line of reasoning is fully consistent with the fact that F 
intraplane coupling does not result in a significant structure relaxation. 

Concerning the electronic states, structural relaxation has serious effects. In contrast to experimental 
findings, previous calculations for the non-optimized geometry resulted in metallic solutions with broad 
conduction bands [19]. This behavior appeared to be characteristic for all spin configurations involving 
parallel spins along the stripes, particularly for our ADFAS, CDFAS, and GDFAS patterns. An insulating 
band gap is observed only for spin configurations involving antiparallel spins along the stripes [18]. While 
the density of states (DOS) for the ADFAS phase is subject to only small changes, the structure optimiza- 
tions for the CDFAS/GDFAS phase yield substantial alterations. The associated spin-polarized partial Ni 
3d DOS, split up into contributions from the four inequivalent Ni atoms, is depicted in Fig. |4] These re- 
sults are based on the LDAh-U scheme, with local interactions on all Ni sites. For the on-site Coulomb 
interaction we set [/ = 8eV and for the Hund's rule coupling Jh ~ 0.8eV [31, 32]. 

Fig. |4] seems to indicate the formation of an insulating band gap, as the consequence of a separation of 
valence bands just above the Fermi energy. Closer inspection of the band structure, however, shows that 
this splitting is incomplete, which could be attributed to the fact that the on-site interaction has not been 
taken into account in the structure optimizations, for technical reasons. Nevertheless, it is clear that the 
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total energy gain (mRyd) 



C/(eV) 


Jh (eV) 


ADFAS 


CDFAS 


GDFAS 








8.7 


19.5 


17.3 


8 


0.8 


163 


156 


157 



Table 4 Total energy gain per Ni atom resulting from the spin-polarization and the simultaneous structural distortion. 
Values are compared for GGA and LDA+U calculations, and for the ADFAS, CDFAS, and GDFAS spin patterns. 

formation of an insulating state is intimately connected to the relaxation of the O octahedra and the stripe 
segregation. 

For the optimized crystal structures, we now compare the spin patterns under consideration from the 
energetic point of view. The energy gain (with respect to the paramagnetic energy calculated for the non- 
optimized La5/3Sri/3Ni04 structure) of the magnetic solutions is summarized in Table |4] In addition to 
the solutions displayed in Fig. HI it seems to be possible to obtain metallic solutions with similar energies. 
However, they are more indicative of a poor convergence than a physically relevant result. All magnetic 
phases are lower in energy than the paramagnetic phase, already for U — 0. In this case, the largest total 
energy gain amounts to 19.5 mRyd per Ni atom for the CDFAS phase. However, on inclusion of the on- 
site Coulomb interaction the energy order is inverted, which contradicts the experimental situation. The 
ADFAS pattern now is favored by 6 mRyd per Ni atom, while the CDFAS and GDFAS configurations are 
almost degenerate. Since the energy differences of the various phases are rather small, this fact again is 
explained by the interplay between the magnetic ordering and details of the crystal structure. Probably, a 
structure optimization within the LDAh-U scheme is required to restore the original energy order 

4 Conclusions 

In conclusion, we have presented DFT calculations for the nickelate La5/3Sri/3Ni04. Optimization of the 
crystal structure results in two characteristic types of lattice distortions: (1) Distinct deformations of the 
NiOg octahedra pave the way for the opening of an insulating gap, which, so far, could not be reproduced 
assuming parallel spins along the stripes. (2) The strong spin-lattice interaction induces a segregation of 
neighboring diagonal filled stripes, as the deformation cooperates with the AF interaction perpendicular to 
the stripes. This new mechanism of stripe segregation is expected to be relevant for a large class of striped 
materials with related crystal structures. Moreover, neither the distortions of the O nor of the Ni sublattice 
further reduce the symmetry initially imposed by the charge pattern. Stripe ordered La5/3Sri/3Ni04 thus 
maintains a Cmmm lattice symmetry, which to verify is an important task for future experiments. Besides, 
there is a clear need for structure optimization on the LDAh-U level, to enable a more reliable discussion 
of the stability of the various phases. Work along these lines is in progress. 
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